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It is shown how on-shell 0(a) improvement can be implemented non-perturbatively in lattice QCD with Wilson 
quarks. Improvement conditions are obtained by requiring the PCAC relation to hold exactly in certain matrix 
elements. These are derived from the QCD Schrodinger functional which enables us to simulate directly at 
vanishing quark masses. In the quenched approximation and for bare couplings in the range < (70 < 1, we 
determine the improved action, the improved axial current, the additive renormalization of the quark mass and 
the isospin current normalization constants Za and Zy. 



1. INTRODUCTION 

The leading cutoff effects in lattice QCD with 
Wilson quarks [|l[ are proportional to the lattice 
spacing a and can be rather large for typical val- 
ues of the Monte Carlo (MC) simulation param- 
eters 1^. Decreasing the lattice spacing at con- 
stant physical length scales means larger lattices 
and therefore rapidly increasing costs of the MC 
simulations. 

An alternative method to reduce cutoff effects 
in lattice field theories is due to Symanzik 
He has shown that the approach of Green's func- 
tions to their continuum limit can be accelerated 
by using an improved action and improved (com- 
posite) fields. A considerable simplification is 
achieved if the improved continuum approach is 
only required for on-shell quantities such as parti- 
cle masses and matrix elements of improved fields 
between physical states [5-7] . 

Symanzik improvement may be viewed as an 
extension of the renormalization programme to 
the level of irrelevant operators. While the 



structure of the possible counterterms is dictated 
by the symmetries, their coefficients have to be 
fixed by appropriate improvement conditions. Al- 
though they can be estimated in perturbation 
theory, a non-perturbative determination of the 
improvement coefficients through MC simulations 
is clearly preferable. 

To 0(a), we have recently carried out the 
non-perturbative improvement programme in 
quenched lattice QCD, thus leaving residual cut- 
off effects of O(a^) only. The detailed results are 
given elsewhere [7-10]. In this report we empha- 
size the general concepts and give a short account 
of the main results. 

We first review on-shell improvement in the 
framework of Symanzik's local effective theory. 
The set-up in finite space-time volume with 
Schrodinger functional boundary conditions is in- 
troduced in sect. 3. We are then in the position 
to state the improvement conditions and discuss 
their evaluation in MC simulations (sect. 4). Fi- 
nally we present our results for the critical quark 
mass (subsect. 4.2) and the current normalization 
constants in the 0(a) improved theory (sect. 5). 



* Talks presented by S. Sint, R. Sommer and H. Wittig at 
Lattice '96, St. Louis, USA, 4-8 June 1996. 
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2. ON-SHELL IMPROVEMENT 

2.1. Lattice QCD with Wilson quarks 

In this section we consider QCD on an infinitely 
extended lattice with two degenerate light Wilson 
quarks of bare mass niQ The action is the 
sum of the usual Wilson plaquette action and the 
quark action 

5F[f/>>] =a^^^(a;)(i^ + TOo)V'(a;), (1) 

X 

where a denotes the lattice spacing. The Wilson- 
Dirac operator 

3 

^^IE[(^*+^m)7m-«V*V^], (2) 

contains the lattice covariant forward and back- 
ward derivatives, and V^. The last term in 
eq. (H) eliminates the unwanted doubler states but 
also breaks chiral symmetry. As a consequence, 
both additive and multiplicative renormalization 
of the quark mass are necessary, i.e. any renor- 
malized quark mass mpt is of the form 

Tor = Z^rriq, rriq = mo - mc, (3) 

where rric is the so-called critical quark mass. 

Chiral symmetry violation is more directly seen 
by studying the conservation of the isovector axial 
current A" . The current and the associated axial 
density on the lattice are defined through 

A^ix) = ^(a;)7^75y^(x), (4) 

P'^ix) = ^(xh^Y^ix), (5) 

where are Pauli matrices acting on the flavour 
index of the quark field. The PCAC relation 

d^A'^^ix) = 2mP«(:r) + 0(a), (6) 

d, = U^t + d,), (7) 

then includes an error term of order a, which can 
be rather large on the accessible lattices [|| . 

The isospin symmetry remains unbroken on the 
lattice and there exists an associated conserved 
vector current. However, it is often advantageous 
to use the current which is strictly local, 

V^{x)^lp{xh^^i^{x). (8) 



The conservation of this current is then also vi- 
olated by cutoff effects, and a finite renormal- 
ization is required to ensure that the associated 
charge takes half-integral values. 

2.2. Symanzik's local effective theory 

Near the continuum limit the lattice theory can 
be described in terms of a local effective the- 
ory i, 

Scfi = So + aSi + 0^82 + ■ ■ ■ , (9) 

where Sq is the action of the continuum theory, 
defined e.g. on a lattice with spacing e ^ a. The 
terms S'fc, k = 1,2,..., are space-time integrals 
of lagrangians Ck {x) . These are given as general 
linear combinations of local gauge invariant com- 
posite fields which respect the exact symmetries 
of the lattice theory and have canonical dimen- 
sion A + k. We use the convention that explicit 
(non-negative) powers of the quark mass m are 
included in the dimension counting. A possible 
basis of fields for the lagrangian Ci (x) then reads 

O3 = TOtr{^^^,F^,}, (10) 

where F^j/ is the field tensor and cr^^ = f [7^1, 7i^]- 
When considering correlation functions of local 
gauge invariant fields the action is not the only 
source of cutoff effects. If 0(a;) denotes such a lat- 
tice field (e.g. the axial density or the isospin cur- 
rents of subsect. 2.1), one expects the connected 
n-point function 

Gn{xi, . . . , X„) = {Z^r iH^l) ■ ■ ■ 0(a^n))co„ (11) 

to have a well-defined continuum limit, pro- 
vided the renormalization constant is correctly 
tuned and the space-time arguments xi,...,Xn 
are kept at a physical distance from each other. 

In the effective theory the renormalized lattice 
field Z^(f)(x) is represented by an effective field, 

(f'csix) = (j)o{x) + a(j>i{x) + a^(j)2{x) + (12) 

where the (pkix) are linear combinations of com- 
posite, local fields with the appropriate dimension 
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and symmetries. For example, in the case of the 
axial current , (fii is given as a linear combina- 
tion of the terms, 

- i^^T'^^5[D, + D^]i;, (13) 

The convergence of G'„(xi, . . . , Xn) to its contin- 
uum limit can now be studied in the effective the- 
ory, 

Gn{xi,...,Xn) = (</'o(a:^i)---0o(a;„))eon 

-ajd^y (0o(a;i)...0o(a;„)/:i(?/))^on 

n 

+ a^{(j)Q{xi) . . . <j>l{xk) . . . 0o(a;n))con 
k=l 

+ 0(a2), (14) 

where the expectation values on the right hand 
side are to be taken in the continuum theory with 
action Sq. 

2.3. Using the field equations 

For most applications, it is sufficient to com- 
pute on-shcU quantities such as particle masses, 
S-matrix elements and correlation functions at 
space-time arguments which are separated by a 
physical distance. It is then possible to make use 
of the field equations to reduce first the number 
of basis fields in the effective lagrangian Ci and, 
in a second step, also in the 0(a) counterterm 0i 
of the effective composite fields. 

If one uses the field equations in the lagrangian 
Ci under the space-time integral in eq. (|l^, the 
errors made are contact terms that arise when y 
comes close to one of the arguments xi, . . . ,a;„. 
Taking into account the dimensions and symme- 
tries, one easily verifies that these contact terms 
must have the same structure as the insertions of 
ipi in the last term of eq. (p|). Using the field 
equations in Ci therefore just means a redefini- 
tion of the coefficients in the counterterm . 

It turns out that one may eliminate two of the 
terms in eq. (p^). A possible choice is to stay 
with the terms Oi, O3 and O5, which yields the 
effective continuum action for on-shcU quantities 



to order a. Having made this choice one may 
apply the field equations once again to simplify 
the term 0i in the effective field as well. In the 
example of the axial current it is then possible to 
eliminate the term Og in eq. ( [T^ ) . 

2.4. Improved lattice action and fields 

The on-shell 0(a) improved lattice action is 
obtained by adding a counterterm to the unim- 
proved lattice action such that the action 5*1 in 
the effective theory is cancelled in on-shell am- 
plitudes. This can be achieved by adding lattice 
representatives of the terms Oi , O3 and O5 to the 
unimproved lattice lagrangian, with coefficients 
that are functions of the bare coupling, go, only. 
Leaving the discussion of suitable improvement 
conditions to sect. 4, we here note that the fields 
O3 and 05 already appear in the unimproved the- 
ory and thus merely lead to a reparametrization 
of the bare parameters go and mo. In the follow- 
ing, we will not further consider these terms. For 
a discussion of their relevance we refer the reader 
to ref. §. 

We choose the standard discretization i^^^ of 
the field tensor Q and add the improvement term 
to the Wilson-Dirac operator (|2|), 

za 

With properly chosen coefficient Csw(5o) this 
yields the on-shell 0(a) improved lattice action 
which has first been proposed by Sheikholeslami 
and Wohlert §]. 

In perturbation theory, the coefficient Cgw has 
been computed to one- loop order by Wohlert [|ll| . 

His result, Csw = 1-1- ciwffo + O(ffo)j with ciw = 
0.26590(7), has recently been confirmed to three 
significant digits 

The 0(a) improved isospin currents and the ax- 
ial density can be parametrized as follows, 

(Ar); = ZA(l + 6Aamq){A^ + acA5^F'^}, 

(Vn); = Zv(l + byam^){V^ + acyd.T;,}, 

{P^r = Zpil + bpam^)P-. (16) 

Here all the fields on the right hand sides have 
been defined in subsect. 2.1 except the tensor den- 
sity T^^ = iil}afj_i,^T°-'il}. We have included the 
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normalization constants ^a.v.p which have to 
be fixed by appropriate normaUzation conditions 
(cf. sect. 5). Again, the improvement coefficients 
&A,v,p and ca,v are functions of go only. At tree 
level of perturbation theory we have = bp = 
by — 1 and ca ~ cy ~ For ca also the 

one- loop result, 



CA = -0.00756(1) X gl + 0{gt>), 
has recently been obtained S . 



(17) 



2.5. The PCAC relation 

We assume for the moment that on-shell 0(a) 
improvement has been fully implemented, i.e. the 
improvement coefficients are assigned their cor- 
rect values. If O denotes a renormalized on-shell 
0(a) improved field localized in a region not con- 
taining X, we thus expect that the PCAC relation 



(18) 



holds up to corrections of order . At this point 
we note that the field O need not be improved for 
this statement to be true. To see this we use again 
Symanzik's local effective theory and denote the 
0(a) correction term in Ooff by 0i. Eq. ( |l8|) then 
receives an order a contribution 



a{{d^{An);{x) ~ 2mn{PRnx)} 



(19) 



which is to be evaluated in the continuum theory. 
The PCAC relation holds exactly in this limit and 
the extra term ( |l9| ) thus vanishes. 

3. THE SCHRODINGER FUNCTIONAL 

3.1. Definitions 

The space-time lattice is now taken to be a dis- 
cretized hyper-cylinder of length T and circum- 
ference L. In the spatial directions the quantum 
fields are L-periodic, whereas in the Euclidean 
time direction inhomogeneous Dirichlet bound- 
ary conditions are imposed as follows. The spa- 
tial components of the gauge field are required to 
satisfy 



U{x,k)\^^^Q = exp(aCfc), Ck = i(f>/L, 



(20) 



with (/) — diag((/)i, 02, 03), and an analogous 
boundary condition with C is imposed at xq = T. 



With the projectors P± ~ ■i(l±7o), the bound- 
ary conditions for the quark and antiquark fields 
read 

P+Mx„=0 = P, P-Mxo=T = P', (21) 

^P_U„^o = P, ll>P+\^,=T = p' . (22) 
The functional integral in this situation |13,nffl. 



Z[C',-p',p'\C,p,p] = [ 

J fields 



(23) 



is known as the QCD Schrodinger func- 
tional (SF). Concerning the (unimproved) action 
S, we note that its gauge field part has the same 
form as in infinite volume. The quark action is 
again given by eq. (Q), provided one formally ex- 
tends the quark and antiquark fields to Euclidean 
times xq < and xq > T by "padding" with 
zeros 0. However, we will use a slightly more 
general covariant lattice derivative, 

Wf,ij(x) = - [X^,U{x, p)iP{x + afi) - i^ix)] , (24) 
a 

where Ao = 1 and — eKp{ia6/L). This modifi- 
cation of the covariant derivative is equivalent to 
demanding spatial periodicity of the quark fields 
up to the phase exp(i6'). We thus have the an- 
gle 9 as an additional parameter that plays a role 
in the improvement condition for the coefficient 
ca|]. 

We are now prepared to define the expectation 
values of any product O of fields by 



(0) = |1/ Oe- 

L ^ J fields 



(25) 



p'=p'=p=p=0 



Apart from the gauge field and the quark and 
anti-quark fields integrated over, O may involve 
the "boundary fields" at time xq = 0, 



C(x) 



<5p(x)' 



C(x) 



<5p(x) 



(26) 



and similarly the fields at xq = T. Note that 
the functional derivatives only act on the Boltz- 
mann factor, because the functional measure is 
independent of the boundary values of the fields. 
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3.2. Continuum limit and improvement of 
the Schrodinger functional 

and 



Based on tlie work; of Symanzik |13 



explicit calculations to one-loop order of pertur- 
bation theory one expects that the SF 
is renormalized if the coupling constant and the 
quark masses are renormalized in the usual way 
and the quark boundary fields are scaled with 
a logarithmically divergent renormalization con- 
stant. 

As in the case of the infinite volume theory 
discussed in subsect. 2.2, the cutoff dependence 
of the SF may be described by a local effective 
theory. An important difference is that the 0(a) 
effective action Si now includes a few terms local- 
ized at the space-time boundaries ||^ . Such terms 
then also appear in the 0(a) improved lattice ac- 
tion. However, by an argument similar to the 
one given at the end of subsect. 2.5, it can be 
shown that they only contribute at order a? to 
the PCAC relation and the chiral Ward identity 
considered in sect. 5. In the calculations reported 
below, the inclusion of boundary counterterms is, 
therefore, not required. 

4. IMPROVEMENT CONDITIONS 

A non-perturbative determination of Cgw and 
ca starts from considering two different lattice 
artefacts, i.e. combinations of observables that are 
known to vanish in the continuum limit of the 
theory. One then imposes the improvement con- 
ditions that these lattice artefacts vanish for any 
finite value of the lattice spacing. The specific val- 
ues of and ca, where this is achieved then de- 
fine the 0(a) improved theory for that value of a. 

Of course, there is a large freedom to choose 
improvement conditions and - on the non- 
perturbative level - the resulting values of Csw 
and CA depend on the choices made. The corre- 
sponding variation of Csw, ca is of order a. This 
variation changes the effects of order a?' in phys- 
ical observables computed after improvement. In 
principle this is irrelevant at the level of 0(a) im- 
provement. Nevertheless, one ought to choose 
improvement conditions where such terms have 
small coefficients. The improvement conditions 
derived from the SF can be studied in pertur- 



bation theory. Such a study provided essential 
criteria for our detailed choice of improvement 
conditions. A further criterion is the statistical 
accuracy that can be obtained in the MC compu- 
tations of the lattice artefacts. 

4.1. Determination of Csw and ca 

Using the operator 



0"^a6^C(y)75^C(z), 



we define the bare correlation functions 
/a(xo) = -\{A^^{x)0^), 
/p(.To) = -\{P\x)0'^). 



(27) 



(28) 



The same correlation functions, but with the 
boundary values C and C" interchanged, will be 
referred to as 5^ and ^p. 

In terms of /a and /p the PCAC relation for 
the unrenormalized improved axial current and 
density may be written in the form 



2/p(xo) 



(29) 



We take this as the definition of the bare current 
quark mass m. The renormalizcd quark mass mp 
appearing in eq. (|l^) is then given by 



mp = m 



Zp{l + bparriq) 



(30) 



At fixed bare parameters, mp and hence also the 
unrenormalized mass m should be independent 
of the kinematical parameters such as T, L and 
a;o. This will be true up to corrections of order 
a^, provided Csw and ca have been assigned their 
proper values. The coefficients may, therefore, be 
fixed by imposing the condition that m has ex- 
actly the same numerical value for three different 
choices of the kinematical parameters. 

In detail we set = 0, T = 2L, {(j)^, (1)2, (j)^) = 
i(-7r,0,7r) and (<^'i, 0^, (/)(,) = i(-57r, 2^, Btt). 
Using the shorthand notation m = rf{xo) + 
CASf{xo) for eq. (|29|), and, correspondingly, m = 
rg{xo) + CASg{xo) for the same equation with gj^ 
and gp, our improvement conditions are 

r/(xo) + CASfixa) = rg{xa) + CASg{xo), (31) 
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Figure 1. Quark mass dependence of Mi for f3 = 
6.4 and Csw = 1.776625. The position of m = 
is marked by the dashed hne. The hopping 
parameter is defined as k = 1/(8 + 2amo). 



at both xo = T/4 and a;o = 3T/4. 

In order to determine Csw we ehminate ca in 
eq. (|3l] ) and consider the combination 

m/= r/(|r)~r,(|T)+7[sK|r)-.,(|r)],(32) 
where 

7 = [r-/(ir) - r,(iT)]/[s,(ir) - Sfi\T)]. (33) 

The equation Mi ~ may then be used to de- 
termine Csw for each value of go. In practice 
we suppress the 0(a) uncertainties in Csw by a 
power of by equating Mi to its tree level value, 
M^^-'la/i, in the improved theory (at the finite 
value of a/L). We further fixed a/L = 1/8, where 
both a good statistical accuracy was achieved and 
the 0(a) uncertainties of Csw are only around 2% 
in a perturbative estimate. 

In general we evaluated Mi for vanishing mass, 

m^rf{T/2) + jsj{T/2). 

As demonstrated in fig. |l|, finding the value of Mi 
at m = is an easy task, since simulations of 
the SF are possible for both positive and negative 
values of the quark mass Q . At each of our nine 
values of go, this is done for at least three values 
of Csw and Mi = M^^^I^/l is solved for Csw by 
linear interpolation of Mi in Csw (cf. fig-|^)- 



In the range < g2 < 1, the results for Csw are 
well represented by 



1 - 0.6565^ - 0.15254 - 0.054.g6 
1 - 0.92252 • 



(34) 



In fig.^ we compare our results to one- loop bare 
perturbation theory and also to tadpole improved 
perturbation theory , for which we have used 



1 + (cW - 1/4)5 



;;2 



(35) 



where 5^ = 59 /uq |Q . Here Uq is the average 
plaquette in infinite volume. 

In a similar way we obtained (0 < 5o < 1) 



ca(5o) = -0.00756 x g^ 



1 " 0.7485g 
1-0.97750'- 



(36) 



Note that the correction term to the axial current 
is rather small and both eq. ( p4| ) and eq. (|3^) devi- 
ate substantially from the one-loop results except 
for values of 59 as small as 5q < 1/2. 

4.2. The chiral region 

As chiral symmetry is violated by lattice arte- 
facts, there is no precise definition of a chiral 
point K = Kc = 1/(8 + 2aTOc) for any finite value 
of the lattice spacing 0| . Rather Kc has an intrin- 
sic uncertainty which is reduced from 0{a^) to 
0(a'^) by non-perturbative improvement. 

Here, we define the value of k, where 

TO, eq. dH), vanishes for C = C = 61 = 0, T = 
2L,xa = T/2. To study the 0{a^) effects we can 
then still vary the resolution a/L. 

Close to Kc the mass to is a linear function 
of K to a high degree of accuracy. The point 



0.004 

„ 0.002 
s 

-0.002 
-0.004 



1.2 



1.4 



1.6 



1.8 



Figure 2. Determination of Csw aX (3 — 6.4. The 
dashed line denotes the tree level value M^'^ \a/L- 
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Figure 3. Non-perturbative improvement coef- 
ficient Csw Results from one- loop bare and tad- 
pole improved perturbation theory are denoted 
by dotted lines and crosses, respectively. 



0.139 



0.138 



0.137 



0.136 




Figure 4. Kc determined for L/a = 8 (circles), 
and L/a = 16 (triangles). The dotted line is the 
one-loop result [|ll],|2|j|] . Crosses represent tad- 
pole improved perturbation theory. 



Kc is therefore easily found by linear interpola- 
tion (or, for /3 < 6.2, slight extrapolation). For 
/3 > 6.8, the results (fig.|) for a/L =1/8 and 
a/L = 1/16 agree on the level of their statistical 
accuracy which is better than 10~^. For lower f3 
small but significant 0{a^) effects are seen. 

Note that the pion mass (in infinite volume) 
will not vanish exactly at k — Kc since there is 
no exact chiral symmetry for finite values of a. 
However, for defined here, the pion 

mass is expected to vanish up to a small lattice 
artefact of order O(a^). 

5. CURRENT NORMALIZATION 

5.1. Chiral Ward identities 

For zero quark masses chiral symmetry is ex- 
pected to become exact in the continuum limit. It 
has therefore been proposed to fix the normaliza- 
tion constants and Zy of the isospin currents 
by imposing the continuum chiral Ward identities 
also at finite values of the cutoff [19-21]. 



In the case of the axial current the relevant 
Ward identity can be written in the form 

da, {x)e'^'^ {A1 {x)Al (y) Q^) = 2i{V^{y) Q^) (37) 

dR 

where the integral is taken over the boundary of 
the space-time region R containing the point y 
and is a source located outside R. In view 
of on-shell improvement it is important to note 



that all space-time arguments in eq. (37) are at 
non-zero distance from one another. 

For the source field in eq. (^ we choose Q°- — 

defined similarly using the primed fields. 

We define the following correlation functions, 
using the unrenormalized improved currents 
[A.Y^ = A- + acA~d,P- and {V,)^ 
acyd^Tf^iy, 

,6 



ya 



/aa(2;o,2/o) = 



\ ^ abc cde 

6L6 ^ 



x(0"^(^i)S(x)(Ai)E;(y)0^), (38) 
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/v(^o) = g^5]*6'^''^(0-m)!;(:.)0^), (39) 

X 

/i = -^^O'^O^). (40) 

At m-q — 0, and using the correct values of 
Csw and ca, a lattice version of the chiral Ward 
identity (|^) is 

ZlfAA{xo, yo) = fi + 0{a^), xo>yQ. (41) 

Compared to eq. ( ^7[ ) we have set v — and in- 
cluded an additional summation over y, thus ob- 
taining the isospin charge. In deriving eq. ( ^ ) we 
have used the fact that the action of the latter on 
the chosen matrix elements can be evaluated due 
to the exact isospin symmetry on the lattice. 

Since the isospin symmetry remains unbroken 
for non-zero quark mass, one need not restrict the 
normalization condition of the vector current to 
the case mq = 0. Using similar arguments as in 
the derivation of eq. ( ^l|) one obtains 

Zy{l + &vamq)/v(a;o) = fi + Oia^). (42) 

Note that the improvement coefficient cy is not 
needed here, because the tensor density does not 
contribute to the isospin charge. 

5.2. Numerical evaluation 

Starting from eqs. (p|,|4^) we impose the fol- 
lowing normalization conditions on the axial and 
vector currents at vanishing quark mass, ruq = 0, 

ZIAa{IT,^T) = A, (43) 

ZvflrihT) = /i, (44) 

where we have set C = C — 9 = 0. It is clear 
that any other choice of the parameters will lead 
to the same result up to effects of the order a^. 
We will use these effects as well as those from 
varying L, T and the insertion points xq, yo as 
an estimate of systematic errors on the normal- 
ization constants. 

Statistical errors on Za and Zy were found to 
be much smaller than those on /aaj /v ^'^^ /i' 
due to the strong cancellation of correlations in 
the ratios /i//aa and fi/fy. A detailed anal- 
ysis of systematic errors at (3 — 6.4 revealed 
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Figure 5. Za and Zv plotted against Qq. Re- 
sults from one-loop bare and tadpole improved 
perturbation theory are denoted by dotted lines 
and crosses, respectively. 



that the spread of results obtained from vary- 
ing C, C", 0, L, T and the insertion points was 
mostly well within the statistical errors. The to- 
tal error was estimated to be at the 1% level for 
Zy and Za- 

Fig.|| shows the measured data for Zy and Za 
obtained for L/a — 8, T/a = 18. They are com- 
pared to one-loop perturbation theory, for which 
we used the results of ref. which have been 
confirmed by one of us. The figure shows that 
contact with one-loop lattice perturbation theory 
can be established for 0-5. At larger cou- 

plings there is also very good agreement between 
our data and one-loop tadpole improved pertur- 
bation theory. 

In the case of the vector current one may 
extend the normalization condition to non- 
vanishing quark masses, such that the ra- 
tio /y(T/2)//i yields the combination Zy{l + 
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bYarriq), from which one can extract the improve- 
ment coefficient by ■ A preliminary analysis shows 
that by is about 40% above its tree-level value of 
by = 1 for 6.0 < /3 < 6.4. 

This procedure can, however, not be applied 
to the axial current, since the mass dependence 
of its normalization is obscured by the physical 
mass dependence arising from the axial density 
in the Ward identity. 



6. CONCLUSIONS 

We have been able to implement on-shell 0(a) 
improvement non-perturbatively in (quenched) 
lattice QCD. The improvement coefficients Csw, 
CA, the critical mass TOc and the current nor- 
malization constants Za and Zy have been de- 
termined for bare couplings in the range < 
(7o < 1. In all cases contact with bare pertur- 
bation theory could be made at couplings around 
(7g « 0.25 — 0.5. The convergence of perturbation 
theory appears to be speeded up significantly by 
tadpole improvement as shown by the crosses in 
figures ||-^. However, for values of /? < 6.8, which 
is the range most relevant for present MC com- 
putations, the quality of tadpole improved per- 
turbation theory is rather non-universal. While 
excellent agreement between perturbation theory 
and our results for Za, Zy is found, substantial 
deviations are observed in the case of Cs„ , ca and 

This work is part of the ALPHA collaboration 
research programme. We thank DESY for allo- 
cating computer time on the APE-Quadrics com- 
puters for this project and are grateful to the staff 
of the computer center at DESY-IfH, Zeuthen, 
where all numerical simulations have been per- 
formed. 
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